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CL|, Trajectories are a central concept in our understanding of classical phenomena and 



< 



also in rationalizing quantum mechanical effects. In this work we provide a way to 

determine semiclassical paths, approximations to quantum averages in phase space, 

f) |. directly from classical trajectories. We avoid the need of intermediate steps, like 

q^ ' particular solutions to the Schroedinger equation or numerical integration in phase 

C^ ■ 

^ ■ space by considering the system to be initially in a coherent state and by assuming 

that its early dynamics is governed by the Heller semiclassical approximation. Our 

result is valid for short propagation times only, but gives non-trivial information on 

\^ ' the quantum-classical transition. 

m ■ 

O : I- INTRODUCTION 

m _ 

Even though wave mechanics does not emcompass the concept of trajectory, this clas- 
^ ! sical notion often appears in the quantum formahsm. On a fundamental level, there are 

^ ■ 1 i n 

trajectory-based theories like path-integrals [1| and Bohmian mechanics [2|] that aim to pro- 
vide interpretations and quantitative rules to deal with general quantum systems. Within 
the orthodox structure of wave mechanics we have the Ehrenfest theorem, which establishes 
the phase-space paths followed by {q)(t) and {p){t), providing an important, but deceptive, 
connection with classical mechanics. An illustration of the continuing interest in the concep- 
tual/experimental status of trajectories in quantum mechanics is a recent implementation 
of appropriate combinations of strong and weak measurements to determine average paths 
of single photons in a two-slit interferometer |3| . 



c^ 



The whole of semiclassics is based on 
mechanical quantities in terms of real [^ 



Newtonian paths. 



6| or complex [7|-. 



;he goal being to express quantum 



14| classical trajectories and their 



stability properties. These approximations are extensively used, e.g., in physical chemistry 
and also in linking the distinct structures of classical and quantum mechanics, where a 
critical point is the understanding of the emergence of classical chaos from a linear theory. 
Usually, in order to calculate semiclassical expressions {q)sc(t) and {p)scit), we first need an 
approximate solution to the time dependent Schroedinger equation \ip{t))sc- The purpose of 
this manuscript is to provide a more direct connection between classical, semiclassical and 
quantum paths, and to discuss the fundamental differences between them. For example, how 
quantum non-locality shows up in the semiclassical trajectories (and vanishes as /i — > 0). 

II. PRELIMINARY REMARKS 

Given a particle of mass fi subjected to a potential V, the equations of motion for the 
expectation values of position and momentum in quantum mechanics, provided by the Ehren- 
fest Theorem, lead to 

where (...)^ = {ip{t)\...\ilj{t)), lipit)) being an arbitrary solution of the time- dependent 
Schroedinger equation. To simplify the notation we suppress the subscript ip hereafter. The 
formal similarity between the previous relation and Newton's second law is quite misleading. 
This is partly due to the widely known fact that 



dV 
dq 
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The two quantities in the above relation tend to coincide only in the limiting case of spatially 
localized wave functions. Some authors refer to the fulfillment of this last condition along 



with Eq. ([T]) as the Ehrenfest Theorem 15|. There is, however, another important difference. 
While in classical mechanics, fix = — dV^/dx is a second order differential equation, relation 
([1]) is an identity, not an equation. Once the quantum dynamics is resolved, \ip{t)) can be 
used to calculate the quantities in ([1]) and to show that both sides of the equality invariably 
coincide. As a consequence, if one resorts to any approximate method to get \4'(t))approx ~ 
\ip{t)), relation ([T]) ceases to hold exactly. In this case the question arises, which side of (P) 
better describes the exact quantum evolution of the mean values ? In principle we are free 



to pick either 

^(gW.o. or --^^^ (3) 

as our trial approximate acceleration (or a linear combination of them). 



III. SEMICLASSICAL EHRENFEST DYNAMICS 

A. The ingredients 

In the present manuscript we will be concerned with semiclassical approximations asc{t)-, 
Psc{t), and qscit) to the corresponding exact quantum expected values. Therefore, our basic 
ingredients are the classical trajectory [qdt) and Pc{t)] and its stability properties given by 
the tangent matrix m. Consider the initial conditions {qo,Po) leading, after a short time, to 
the phase-space point {qf,Pf)- If we take the neighboring initial point (go + ^qo^Po + ^Po), 
one ends up at (g/ + Sqf,pf + Spf). In the linear approximation m is defined by 





iTLpq iTipp 

where h and c are constant parameters, that will be conveniently set, with dimensions of 
position and momentum, respectively. 

In the semiclassical domain it is natural to consider localized states, which uniformly 
tend to a classical point in phase space as /i — ?> 0. For this reason, a canonical coherent state 
is taken as our initial ket, the corresponding wave function being given by the minimum 
uncertainty Gaussian packet 

^o(x) = {x\zo) = ^^T7i^ e ^5^+.P"(- 2 ) , (4) 

where \zo) is a coherent state with zq = (qo/b + ipo/c)/\/2. The position and momentum 
uncertainties are Ag = b/\/2 and Ap = c/v^, with be = h. These relations also provide 
convenient scales b and c used in the definition of the stability matrix. In order to satisfy 
the minimum uncertainty relation we must have b ~ h'^^'^ and c ^ W for any 7. However, 
for 7 7^ 1/2, either b/c or c/b diverges as /i — )■ 0, leading to a non-physical infinite squeezing 
in the classical limit. Thus, we assume from now on that 

b ~ yh and c ~ vh , (5) 



so that the classical limit corresponds to a rescaling, the geometrical nature of the state 
remaining unchanged. 

Finally, given the classical path and initial state, we choose the simplest semiclassical 
approximation to describe the time evolution of \zq), namely. Heller's "Thawed Gaussian" 

dfriqq + IfTiqp 

where Qc and Pc are the classical position and momentum as functions of time, S is the 
classical action, and C = ij^pp — ^i^pq) / {''^qq + i^^qp)- Due to its extreme simplicity, l^^cP 
is a Gaussian for all times and the Heller wave function is usually a poor approximation to 
describe the quantum wave function for longer propagation times. Therefore, our results 
will be valid only on a short-time scale. This limitation, however, will not prevent us to get 
useful information. 



B. Semiclassical trajectories 

Since the above formula is a Gaussian for all times, its use together with the first option 

in ([3]) leads simply to asc = dc- In this case, no vestige of the quantum behavior is left. On 

the contrary, if one employs (jHD and the second option in Q, the result is non-trivial. Thus, 

we define the semiclassical acceleration as 

1 /dV\ 1 /""^ , , , . ModV" 

asc = (-7^/ = / dx\tpsc{x,t)\ -— , (7) 

jj, \aq I ^^ jji J-oo dx 

where ipsc{.x) = {x\i/j(t))sc- This definition provides corrections to the bare classical trajec- 
tory because the Gaussian extends over the space "sensing" the different values taken by 
dy/dx, i. e., the whole force field. It is important to note that the potential cannot be 
singular for any finite x, since the Gaussian is non-zero everywhere. We begin our analysis 
of ([7]) in a very general way, by assuming that the potential is an arbitrary polynomial: 

N 
Vix) = J2 ^nX^ ■ (8) 

In principle, N is supposed to be finite, but we may take A^ — ;■ oo in cases where convergence 
of age can be guaranteed. Substituting the above relation in the general prescription ([7]) we 
get 

asc = > r- / X / dxx e ■'t*' , (9) 



with cr(t) = hJm^ + m^ This leads to 



N 



/^ n=l 



a{t) 



n-\ 



H 



ji— 1 



a{t) 



(10) 



2z 

where H^ denotes the Hermite polynomial of degree k [notice that i'^Kkiw) is a real function] 
and the time dependence was made explicit. By integration we get Pscit) and qsc(t). We 
recall that in order to get these approximations to the quantum expected values one has only 
to solve the classical equations for qdt) and the stability matrix m, not the Schroedinger 
equation. 

C. General properties 

As any consistent semiclassical approximation, age must coincide with the classical and 
exact quantum values for Gaussian packets in linear and harmonic potentials {V = const, 
V (X X, and V oc x^). This is indeed the case. For V = const, formula (^^ gives immediately 
ttsc = 0, while for V = aix, relation (ITOl) gives the constant acceleration age = —ai/fi- 
Finally, for V = a2X^, we get a^c = —2a2qc(t)/fi. Setting ^2 = /iCo'^/2 we obtain agcit) = 
—uj'^qcit), as expected. For higher order potentials the classical, semiclassical, and quantum 
results no longer coincide. 

Before going into more quantitative examples, we establish the general lower order cor- 
rection to the classical acceleration in the general case of an analytic potential. To do 
that we must realize that the Planck constant appears in ([71) [or ( TTOl) ] implicitly, through 
cr oc 6 oc \/Ti. Replacing 

H„(,) ^ g(2,)"-(-l)'^^^ (11) 

into ( ITU]) and collecting the lower powers of a we get 



ttc - 



This shows that our result is consistent with classical mechanics plus a main correction of 
order h^ and smaller terms (involving higher powers of h and higher derivatives of V). Per- 
haps the most important point about the above relation is that these approximate quantum 
paths cannot come from any effective or modified potential, because they are not Hamil- 
tonian. Differently from what happens in classical mechanics, the instantaneous dynamics 
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of a phase-space point depends upon higher order derivatives of the potential. This is a 
trace of quantum non-locahty. While in the Newtonian mechanics the particle only needs 
to "know" its position [to get V{x)] and the force to which it is subjected (ocdl^/dx), both 
local quantities, in the semiclassical case all derivatives of the potential are needed, which 
amounts to knowing the function V for all values of x, no matter its distance to the particle. 
Any change in the potential in a far region instantly influences the semiclassical particle's 
path. Of course, this non-locality vanishes as h —> 0. 

Note that the correction to the classical acceleration is, in general, non-vanishing even 
for t = 0. This can be understood by realizing that our expression yields, by construction, 
the exact quantum result as t — )■ 0. For sufficiently short times the evolution operator can 
be written as / — ^Ht — ^HH^, so that the expectation value of momentum is, to second 
order in t 

= Po4<(|V'M).p])»-|;({p.^)) +0(a 

where { , } stands for the anticommutator. In the quantum mechanical case (where the 
Ehrenfest Theorem holds) we define the exact acceleration as 

Id,.,,, 1 /dV\ t /i.d^V]\ ^,0, 

a,.„„, = --m) = -- ( g-)^__ - ^ ({p. g|r})^__ + 0(e) , (14) 

which is asc for t = 0. This guarantees a consistent quantum-classical interplay in the 
sense that the semiclassical expression tends to become classical for ^ — )■ while it is the 
exact quantum result for t — > 0. This closes our discussion for short-time propagation in 
ID analytic potentials. In the next section we illustrate the procedure by applying ( TTOl) in 
the case of a cubic potential. In addition, we give an example involving a non-analytical 
potential, for which expression (fTOj) is not valid and ([7]) must be used directly. 



IV. APPLICATIONS 

A. Cubic potential 

Here we consider in more detail the classical, semiclassical, and quantum accelerations in 
the particular case of a cubic potential V = ax"^. The quantum result up to first order in t 
is 



a. 



quant 



3a 






3a 



% + j + -goPo 



Let us calculate age also up to first order in t. The full semiclassical expression reads 



3a 



qdtf + 



a 



(tf 



(15) 



(16) 



which coincides with (TT2|) for a cubic potential. We have qdt) =10+^ ^ + C'(t^), from 
which it is easy to show that mqq ~ 1, mqp ~ t^, so that cr^ = 6^ + O(t^). Substituting these 
expansions into (TT6|) we get 

Thus, in this case, the quantum and semiclassical results coincide not only to zeroth order, 
but also to first order in t. On the other hand we know from (1121) that age goes to Oc as 
/i — 7- 0. We can immediately write 





3a 


\ b^ 2t 


Clsc ~ 





% + ^ + -qm^ 




/i 



(17) 



3a 



2 2t 
% + —QoPo 



(18) 



since 6^ ~ h. In particular, it is interesting to note that the classical equilibrium point qo = 
and po = 0, leading to a^ = 0, does not occur in the semiclassical and quantum cases, since 
there is always a "residual" acceleration age = Equant = 3a6^/2yU. 

B. Wave packet impinging on a step potential 



A crucial process in wave-packet dynamics is the collision with a potential wall. Here we 
examine the evolution of (g) for the Gaussian packet (HI) bouncing off a step potential. 



V{x) 



, for X < 
Vq , for X > 



(19) 



for the cases where po < -\/2/rVo. From the classical trajectory, 



ge(t) = <^ ^ (20) 

2f(t-to), for t>to, 

with to = f^'Qa/poi we extract the two needed tangent matrix elements: 

f 1 , for t < tn 
m,,{t) = ' - ° (21) 

I -1 , for t > to 
and rjigpit) = htmqq{t) / fib'^ . The semiclassical acceleration, as defined in ([7]), is simply 

(22) 





(^sc 


--|^.c(0,t)2 




3]) and (21j), 






asc{t) = 


Vo 
/i0ra(t) 


e V / , 

(p(.t-to)Y 


for t < to 
for t > to , 



(23) 



where cr(t) = b\ i-j^) + 1. Numerical integration leads to Qscit), which is depicted in figure 
1 along with qdt) and the time evolution of the exact quantum average {q){t). The initial 
position and momentum are qo = and po = 1, and the barrier is located at g = 1 with 
height l^ = 5, in a system of arbitrary units. The other parameters are b = 0.1 and fi = 1. 
In the top panel h = 0.05 and in the bottom panel h = 0.1 (Note that the classical action 
for these trajectories is of order of 1). The time interval in the horizontal axis ends shortly 
after the classical turning time. We conclude that even for this discontinuous potential the 
semiclassical results are quite consistent for times shorter than the classical turning time. In 
particular, the tendency of the quantum turning point to recede is seen in the semiclassical 
plots. In fact, this effect may be artificially enhanced in our approximation because in the 
quantum dynamics the packet deforms such that its penetration into the barrier is smaller 
in comparison that of the semiclassical Gaussian packet. The deformation is stronger for 
larger values of h (see the bottom plot). 

V. CONCLUSION AND PERSPECTIVES 

We have developed a direct way to determine semiclassical trajectories for short propa- 
gation times. Despite this limitation we were able to reproduce some important quantum 
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features, like the tendency of turning points to recede in the quantum regime. In addition, 
we identified how quantum non-locahty manifests itself in the semiclassical phase-space, de- 
stroying the original Hamiltonian structure. In this context it is not possible to resort to 



any extended Hamiltonian formalism [1?| (see also[18|) . 

Perhaps, a promising perspective to be drawn from this approach is related to its gen- 
eralization to higher dimensional problems in connection with chaos. On the one hand we 
have classical mechanics which presents sensitivity to initial conditions in a variety of sys- 
tems, and, on the other hand is quantum mechanics, whose linearity prevents exponential 
separation of (r) and (r + dr). In principle, our semiclassical paths should be in the middle- 
way, which raises the question on the behavior of the Liapunov exponents associated to 
semiclassical paths and on their dependence on h. 
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in the right figure h = 0.1. The time interval in the horizontal axis ends soon after the classical 
turning time. 
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